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Abstract 

A non-linear backward equation with diffusive terms is postulated for 
the probability density that depends on the Bohmian quantum potential. 

An associated nonlinear Schrodinger equation is also introduced and ex¬ 
tension of the analysis to several particle compounds is sketched along 
with the implications following from this formalism regarding the non¬ 
conservation of probability in the non-equilibrium regime. Some further 
conclusions are educed with respect to the Onsager’s relations and the 
generalized optical theorem. 

In one of the early papers on wave mechanics, Darwin [1, 2] examined an 
interesting example of a two-particle quantum interference involving a set up of 
a coherent pair of particles. A few years later interparticle quantum correlations 
between two bodies became notable with the work of Einstein, Podolsky and 
Rosen[3, 4] and since then have been at the heart of quantum mechanics. The 
issue of non-separability between interacting particles had been particularly 
emphasized by Einstein in his own version of the EPR thought experiment 
[5]. Schrodinger in a similar vein claimed quite succinctly that according to 
quantum theory maximum knowledge of a composite system does not necessarily 
mean fullest knowledge of the individual particles [6, 7]. This type of reasoning 
that emphasizes the influence of the particular experimental conditions for the 
available type of possible predictions has been extensively used by Bohr in his 
reply to Einstein-Podolsky-Rosen [4] and is best illustrated in the ingenious 
thought experiments devised earlier by him and Rosenfeld for the measurability 
of fields [8]. More recently Mermin in particular emphasized the physical reality 
of these subsystem correlations in comparison to the correlata underlying them. 
From a rather different perspective, the quantum potential which has a pre¬ 
eminent role in Bohmian dynamics, is associated with the interrelatedness in the 
quantum world and suggests a quality of irreducible wholeness in the domain 
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of quantum phenomena [10, 11]. Here we will modify the continuity equation 
for the probability density by including a term proportional to the quantum 
potential and examine some of its consequences. 

Another issue that needs to be addressed is the time asymmetry observed 
in the physical world in contradistinction to the time symmetrical form of the 
fundamental equations of evolution [12]. The resolution of the contradiction in 
classical mechanics is usually related to an asymmetry in initial conditions. In 
quantum mechanics in particular, Sclrrodinger equation becomes time reversible 
if we make the identification of the probability amplitude with the time sym¬ 
metrical complex conjugate one. It is generally appreciated that the quantal 
irreversibility is the result of initial conditions or the outcome of the interaction 
of a few particle system with a reservoir but recently it has been shown experi¬ 
mentally [13] that even a few particles may exhibit signs of irreversibility. The 
collapse of the wavefunction on the other hand is generally accepted to impose 
a genuinely time asymmetric element in the quantum formalism. 

In this work we follow a different approach, by demonstrating that the quan¬ 
tum potential term mentioned above may be employed in a way that implies that 
the evolution equation for the probability density is genuinely time-asymmetric. 
It deserves notice that the equation that will be proposed depends on a final 
value condition so it turns out to be backward [14] and resembles on that respect 
the Kolmogorov backward equation [15], the adjoint transport equation [16] or 
the Black-Scholes one known from mathematical finance [17]. Following this we 
extend our study for the multiparticle and present some conclusions related to 
the quantum equilibrium hypothesis, Onsager’s relations, the generalized optical 
theorem and finally we recapitulate our findings in the final section. 

Generalized continuity equation and its associated non-linear wave equation. 
In a recent paper [18] it has been argued that quantum theory may not be ex¬ 
tended straightforwardly to the macroscopic domain. This could lead us to seek 
generalizations of the present quantum formalism. Since standard non-linear 
versions of the Schrodinger equation [19, 20] seem to face certain difficulties 
[21] we try instead to generalize the Bohmian-de Broglie formulation. Employ¬ 
ing the position representation, we consider initially a neutral, spinless particle 
of mass ?n, assuming that the time varying potential is V(t). The fundamen¬ 
tal quantum mechanical equation of motion for the single particle probability 
amplitude reads as 
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The familiar polar substitution 'I' = Re lS ^ h gives the Hamilton-Jacobi-Madelung 
equation when we equate the real parts 
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We may also define the quantum potential as 
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which is important in Bohm’s version of the de Broglie-Bohm mechanics, al¬ 
though it is possible to formulate the theory without any reference to it. Equat¬ 
ing the imaginary parts we obtain 

1 

— + — (2 VR ■ VS + RV 2 S) = 0. (4) 

at 2m 

Born’s rule which at least in text-book quantum mechanics is considered a sep¬ 
arate axiom of quantum physics provides a conceptual and calculational bridge 
between theory and experimental results. Taking into account then the identi¬ 
fication p = u 2 for the probability density and multiplying both sides with 2 u 
yields the familiar continuity equation 


dp 


VS 
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From (2) and (4) we may determine in principle the action S and the amplitude 
R. We seek to modify the previous expression (4) including a diffusion term and 
a linear term proportional to the quantum potential as follows 


^ + — (2Vu • VS + uV 2 S ) = (—— ) ( V 2 u - . (6) 

This formula may be considered as a non-equilibrium extension of (4) in the 
sense that as u — > R the former reduces to the familiar equation that holds at 
quantum equilibrium. The given expression remains time-asymmetrical which 
immediately implies non-unitary evolution (the inverse proposition is not true 
nevertheless). The negative coefficient suggests a backward evolution for u(t). 
We conclude then that not only the evolution equation is not time-invariant 
but that also there is time asymmetry in the imposition of terminal instead of 
initial conditions. An element then of finality is introduced in (6) since it relies 
on terminal values evolving backwards in time. The homogeneity condition 
is satisfied since for an arbitrary solution u and k a real number it follows 
immediately that ku is also a solution. Hence the rescaling transformation v! —» 
ku constitutes a similarity transformation that leaves the equation invariant. An 
imaginary quantum potential would give rise to dissipative terms that do not 
preserve unitarity in the general case. At the same time we must ensure that 
the expression reduces to the standard form at equilibrium. This suggests an 
extension of (1) written as 




V 



[ v 2 \AW 

{ #F 


J 


(7) 


where $ = ue lS ^ h and the wavefunction d' = Re lS ^ h is governed by (1) with 
appropriate final conditions. The existence of a distinct wave function with 
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the same phase as but different real amplitude is reminiscent of the double 
solution theory put forward by de Broglie [22], but it must be noted that in (7) 
$ propagates backward in time in the configuration space in the same manner 
as so it is not a singularity in the real three dimensional space as de Broglie 
envisaged. Dynamical terms of the form given in (7) have also been used in 
the past in the real but not in the imaginary part of the Hamiltonian [23]. It 
must be underscored that when the terminal condition is fulfilled, which means 
u(t —► oo) = R{t — > oo) so —> oo) = ^(t —> oo) the superposition principle 

holds and we recover (1) from (7) even though the coefficient K/m is obviously 
different from zero. Hence the magnitude of the non-linearity coefficient may 
be non-zero nevertheless there is an appropriate limit that we recover (1). On 
the other hand a nonlinear Schrodinger as in [19] reduces to the linear one only 
when the nonlinearity constant becomes zero. Finally, (7) is time irreversible 
even when $ —> $* and respects the homogeneity condition [20] which says that 
if $ is a solution of (7) so is —> Zd>, where Z a constant complex number. 
Based on (7) we can present extensions of standar quantum formulas. If we 
write for example the expectation value of an operator A as 

(i) = J <S>*A<I>d 3 x, (8) 

its time development is given by 


d(A) 

dt 
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where H the Hamiltonian and of course the first two terms make up the familiar 
expression in the Schrodinger picture whereas the last term reflects dissipative 
effects. According to (6) and taking into account that p = |<I>| 2 , the correspond¬ 
ing expression for the probability density is 


dp 

dt 
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The left part is familiar from the continuity equation but the right non zero part 
implies non-unitary evolution. The existence of these terms suggests a possible 
interplay between spatial-temporal coherence and diffusion phenomena. Diffu¬ 
sion terms in the continuity equation have been also been added in [24, 25] (see 
also [26, 27]). In addition equation (10) remains invariant when R —» kR and 
p -* A p. It must be noted though that the diffusion constant in our proposal 
is negative, namely —h/m as we mentioned. Without the the nonlinear and 
the quantum potential term, (10) is a non-equilibrium backward Fokker-Planck 
type equation that conserves probability locally. The inclusion of the quantum 
potential and nonlinear terms destroy this property and is responsible for its 
non-conservation. It is evident that superluminal communication is precluded 
by the very fact that the relevant equation of motion for the probability density 
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depends on the imposition of a final condition. This preferential status of ter¬ 
minal conditions instead of initial means that we could use (6),(7), or (10) to 
make retrodictions but not predictions. It may be commented that this special 
importance of final conditions can be clearly illustrated with the help of the well- 
known thought delayed choice experiment [28] and even the Bohrian emission of 
photons during atomic transitions in the old quantum theory era. At the equi¬ 
librium limit these reduce to the more familiar Bohmian equations so they may 
been seen as generalizations of the latter in the non-equilibrium domain. The 
presence of the diffusion terms would tend to make u smoother but the quantum 
potential driving term and especially the highly nonlinear midterm could have 
the opposite effect. On physical grounds we would expect irregularities in p(x, t) 
to become less smooth with backward time evolution which may cause problems 
for the behaviour of the solutions like arbitrary rapid growth in finite time [14]. 
The presence then of the nonlinear terms is of special importance, since if they 
tend to create irregularities as we propagate backwards in time and dominate 
over the decay term that will decide in favour of the ill-posedeness of the prob¬ 
lem and would demand sophisticated regularization techniques. Obviously, we 
could also include advanced vector potentials and even noise terms in (10) [29]. 
The latter ones require special care though because of certain mathematical 
subtleties involved [30]. 

It must be emphasized though that there are no non-intersecting particle 
trajectories apart from the case of equilibrium when (8) becomes the familiar 
density current equation and the coefficient sets time scale equal to h/m. We 
now seek a guidance equation modifying the de Broglie-Bohm one. Motivated 
by (11) we write a formula for the velocity v of a single particle that would 
accompany (7) as 
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(ii) 

The first term in the above formula is the velocity term —, the second one 
proportional to VZ?x| <4> | 2 is an osmotic type contribution to velocity according 
to [10, 31] (note the negative sign of the diffusion coefficient) and the third inte¬ 
gral term expresses the non-local character of this generalized guidance equation 
and clearly violates the uniqueness propery of trajectories. Since the last two 
term cancel each other when $ = 4/ we recover the standard guidance equa¬ 
tion. Having established the basic elements of the formalism (relevant evolution 
equations) we proceed by examining some further developments and examples. 

It is interesting to note, as it is clearly seen in (10) or (11), that even when 
the gradient of the action is zero, the probability density may be nonstationary 
which implies the possibility of backflow effects. The inverse is also true, (10) 
allows a steady state density even when the the momentum is non-zero. Since 
ft/m is very small we expect that these deviations from (5) to be very weak. 

Seeking a similarity solution we proceed by defining the scaling anzatz rj = —j= 
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and obtain a nonlinear ordinary differential equation which still may be not 
amenable to analytic solution. Instead we could make the substitution p = e 2h 
to (10) when the spatial derivative of the phase is considered negligible and find 


dh 

dt 



V 2 h + (V/i) 2 + ~^Q 


( 12 ) 


The above resembles the celebrated Kardar-Parisi-Zhang equation [32] that has 
many applications in the field of surface growth. The crucial difference is that 
it propagates backward in time and the last quantum potential driving term 
should be considered deterministic in general instead of stochastic. 

Non-conservation of total probability and entropy production rate. The rate 
of change of the total probability is expressed as 


dP d f d 3 xp f 3 dp / h 

dt dt J X dt \ m 
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The second term makes apparent a connection to Fisher information, which is 
a measure of system disorder [33, 34]. As p aprroaches R the rate of change of 
the total probability becomes zero. Making use of the identities = \7hiR 
and = V 2 lnR + ( VlnR ) 2 and allowing p —> R 2 + 5p , the variation of the 
total probability rate SI may be written as 


Differentiating over time the Gibbs entropy G = — f d 3 xp\np which relies on 
the spatial distribution of the probability density gives the following rate for the 
production of entropy 


^V 2 - V/mRV - S7 2 lnR) 6p. 


(14) 
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, (15) 

The second integral remains unaffected when time is reversed and can be iden¬ 
tified as a measure of entropy production S and a sign of system irreversibility. 
The other term in (15) is associated with entropy flow that follows from (5). 
Again near the linear approximation regime p —> R 2 +6p the variation of entropy 
production is given by 


dG 

dt 


J d 3 x + — J d 3 x(l + Inp) 


5Y, = VFiUjVFj. (16) 

i,j 

In the last equality, F l :] are generalized thermodynamical forces near equilibrium 
[35] and Lij the corresponding Onsager coefficients. Since microscopic time 
reversibility and the resulting detailed balance do not hold it follows according 
to (15),(16) that these coefficients may not be symmetrical. 
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Quantum scattering and the generalized optical theorem. Up to this point we 
have only considered time dependent evolution equations. It is possible never¬ 
theless to make use of a time-independent version of (7) invoking its imaginary 
potential. The time-independent variant of this nonlinear wave equation for a 
particle with energy E, mass m and wave-vector k is written as 


(V 2 + k 2 - U r ) 


( v 2 ^Wv F 

V vF$MF 


VWFTF ) 


<F(r) = 0 


(17) 


where k =, U{ r) = ^pU(r) and U(r) a central potential. When <f>(r) = 'l'(r) 
the formal solution of (17) can be written as 


'F(r) = 'Fg(r) + J G k (r, r')U(r')^/(r')d 3 x (18) 

where G k (r, r') the corresponding Green function and ’Fo(r) a solution of the 
homogeneous equation (V 2 + fc 2 )To(r) = 0 [36]. Due to the nonlinear character 
of (17) it is not possible to write a relation similar to (18) for < f ) (r). We can 
proceed though by assuming that at large distances from the potential region 
the solution <h(r) has the asymptotic form 

$ k ^e ikr +— /(k,k'). (19) 

r 

Following the analysis by Glauber [37, 36] for two scattering vectors k, k' (17) 
produces the identities 


$-^V 2 $ k - $ k ^V 2 $_ k = 0 (20) 

and 
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Integrating (20) and using Green’s theorem over the volume of a very large 
sphere gives /(k'.k) = /(—k', —k) which expresses the reversibility of the scat¬ 
tering processes/amplitudes for two directions [37]. Furthermore, if the scatter¬ 
ing potential has inversion symmetry then the scattering amplitude is symmetric 
/(k',k) = f(k, k') which is deary related to the principle of detailed balance 
[37]. This conclusion nevertheless does not contradict our previous discussion 
which relied on a time-dependent evolution equation. Following the same pro¬ 
cedure for (21) we find 


±[f(k',k)-r(k,k')} = ^ J /*(k",k)/(k",k)div 


i 

47T 


vFW 


v 2 #iF \ 
vW ) 


$ k $£,d 3 ;r. 


( 22 ) 


7 



Setting k = k 7 gives 


Im/(k, k') 
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The integrand in the right hand of (22) is easily identified with the divergence of 
the current density j = |4>k| 2 -yf- Writing <I>k = e 9 ^ r ^4'k (since they both have 
the same phase) where g(r) an unknown real function to be determined yields 
the integral equation 


47rlm/(k, k') = ka 


V 2 5 (r) + V<?(r) • VZn^FF + (V</(r)) 2 


e 2 9( r ) |vj/ k | 2 d 3 


(24) 

It is clear that we recover the optical theorem for elastic scattering when $k = 
Tk which it is known that expresses the particle conservation requirement but 
more generally deviations may be observed according to (23) or (24). When 
the above quantity is positive attenuation of the incident beam is implied. It 
is interesting to note that the presence of the imaginary potential in (17) intro¬ 
duces non-unitarity since the total particle number is not conserved, but does 
not affect the reversibility of the scattering amplitude [37]. 

Many particle case and conclusions. It is straightforward to extend our 
previous results for the case of a compound of several interacting particles of 
equal masses m 



The presence of the many particle quantum potential in the continuity equa¬ 
tion signifies an important feature, namely the dependence of the multiparticle 
density on a quantity that exhibits non-local characteristics (but since (25) is 
backward in time there is no possibility of non-local signalling). The intercon¬ 
nectedness of the system constituents is expressed then through the form of this 
non-local quantum potential. The constituent particles then are reidentified 
and incorporated into a substantial whole unified by a single unifying principle 
and cannot been considered as distinct aggregated components. This means 
that the component parts even though are retrievable in principle, they are not 
intrinsically unaffected by their union [38, 39, 40]. Similar considerations in a 
broader sense have also been expressed various authors. Whitehead for example 
has spoken about the concept of a dynamic and contextual organic whole [41], 
as of course Bohm and Hiley [42, 43] already mentioned. 

In summary, in the present genuinely multiparticle formulation ontological 
priority is given not to the monadic particle but to the emergent several-particle 
composite. We conclude that the presence of the many particle quantum po¬ 
tential in the modified backward continuity equation suggests non-separability 



and illustrates the irreducible nature of a multi-particle compound. It is rea¬ 
sonable to deduce then that we cannot examine each particle singulatim neither 
we should consider the quantal multiparticle as a mere mereological aggregate 
of distinct co-present entities accidentally unified, but rather as a substantial 
unified whole that constitutes the production of a new physical reality. This 
particular formalism also ensures that the time anisotropy inferred from the 
evolution equation and its dependence on terminal values is nomological and 
intrinsic in nature rather than contingent and accidental. 
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